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Abstract: We investigate the dynamics of a neutral and a charged particle around a black 
hole in modified gravity immersed in magnetic field. Our focus is on the scalar-tensor-vector 
theory as modified gravity. We are interested to explore the conditions on the energy of the 
particle under which it can escape to infinity after collision with another neutral particle 
in the vicinity of the black hole. We calculate escape velocity of particle orbiting in the 
innermost stable circular orbit (ISCO) after the collision. We study the effects of modified 
gravity on the dynamics of particles. Further we discuss how the presence of magnetic field 
in the vicinity of black hole, effects the motion of the orbiting particle. We show that the 
stability of ISCO increases due to presence of magnetic field. It is observed that a particle 
can go arbitrary close to the black hole due to presence of magnetic field. Furthermore ISCO 
for black hole is more stable as compared with Schwarzschild black hole. We also discuss the 
Lyapunov exponent and the effective force acting on the particle in the presence of magnetic 
field. 
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I. INTRODUCTION 


Theories of modified gravity (such as f(R) theory, Lovelock gravity, Gauss-Bonnet theory etc) 
are constructed by adding curvature correction terms in the usual Einstein-Hilbert action through 
which the cosmic accelerated expansion might be explained [l| (see also jij] for reviews on modified 
gravity). Such correction terms give rise to solutions of the field equations without invoking the 
concept of dark energy. To find the dynamical equations one can vary the action according to the 
metric. There is no restriction on the gravitational Lagrangian to be a linear function of Ricci 
scalar R [3]. Recently some authors have taken into serious consideration the Lagrangians that 
are “stochastic” functions with the requirement that it should be local gauge invariant [4], This 
mechanism was adopted in order to treat the quantization on curved spacetime. The result was 
that corrective term in the Einstein Hilbert Lagrangian arises due to either background geometry 
and interactions among quantum fields or gravitational self interaction [5]. Furthermore, it is also 
realized that such corrective terms should be incorporated if one wants to obtain the effective 
action of quantum gravity on Planck scale Q]. Besides fundamental physical motivation, these 
theories have acquired a huge interest in cosmology as they exhibit inflationary behaviors and that 
the corresponding cosmological model seem very realistic [7, 8]. In this article, our focus will be on 
the scalar-tensor-vector theory (will be referred as MOG) and the Schwarschild-MOG black hole 
(MOG) 



the dynamics of a particle (either massive or massless) around the gravitational source such as a 
black hole (BH) is important because it is responsible for understanding the geometrical structure 
of spacetime near the BH. Geodesics may display a rich structure and convey a very reliable 
information to understand the geometry of the BH. There are many types of geodesic motion 
but the circular geodesics are specially important. The exponential fade-out of a collapsing star’s 
luminosity can be explained by the circular geodesics as given in 12, [13[]. The motion of test 
particles helps to study the gravitational fields of objects experimentally and to compare the 
observations with the predictions about observable effects (light deflection, gravitational time delay 
and perihelion shift). 

In the surrounding of the BH, a magnetic field is generally present [14|], due to the presence of 
plasma in the vicinity of the BH. The accretion disk or a charged gas cloud is primarily responsible 
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for the magnetic field field Jl5, 16]. The magnetic field is stronger in the vicinity of BH’s event 


horizon however it does not effect t 


moving around a BH is effected 


13 , 


re geometry of the BH, but the motion of the charged particle 
18]. The magnetic coupling process is likely responsible for the 


stability of black hole with its accretion disc 


191 ]. According to this process, angular momentum 


and energy are transferred from the black hole to its surrounding disc. The magnetic field plays an 
important role in transferring sufficient energy to the surrounding particles for escaping to spacial 
infinity [20, 21]. Other interesting processes around BHs include evaporation and phantom energy 


accretion onto BHs 


22], In this article, we revisit the model of Zahrani et al [23|] for a MOG black 


hole and explore the effects of modified gravity. It involves the collision of a bounded particle 
with an unbounded particle in the vicinity of the black hole. The main interest lies in finding the 
conditions of escape of a bounded particle after the collision. 

The outline of the paper is as follows: In section II we develop the basic equations and then 
derive an expression for escape velocity of a neutral particle. In section III we discuss the strength 
of magnetic field, and the equations of motion of the charged particle moving around weakly 
magnetized MOG and escape velocity for particle also calculated. Force acting on the particle is 
studied in section IV and geodesics of the particles moving around the MOG are discuss in section 
V. Lyapunov exponent is explained in section VI. In section VII and VIII, trajectories for effective 
potential and escape velocity of the particle are presented respectively. Conclusion is given in last 
section. We will study the motion of particle in the equatorial plane to simplify the calculations. 
Throughout this work we use the following metric signature and assume c = 1. 


II. NEUTRAL PARTICLE DYNAMICS AROUND MOG BLACK HOLE 

Motion of particles around a central massive object under the affect of a central force is a 
well-studied problem of classical mechanics (or rather Newtonian mechanics). In particular we can 
think of the following problem in the present context: consider a particle of mass m moving in a 
circular orbit around another object of mass M such that M » m. Now for the particle to escape 
from the gravitational field of M, particle’s initial velocity must be more than the escape velocity. 
The particle can gain the escape velocity either from an external force acting on it or by hitting 
(or colliding) a test particle with the particle in circular orbit. Since the collision leads to transfer 
of energy as well, the particle will escape from the circular orbit if its energy after collision is more 
than a critical energy or escape energy. If however, the energy of the particle after collision is small 
than the critical energy, the particle falls towards M. 
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The relativistic version of the above scenario was investigated by Zahrani et al 


23j. They studied 


the motion of a charged particle in the vicinity of a weakly magnetized Schwarzschild black hole 
and focused on the bounded trajectory lying in the black hole equatorial plane. For the charged 
particle to escape from the innermost circular orbit, another particle (which is neutral and coming 
from sufficiently far distance) hits the charged particle. The authors obtained the corresponding 
conditions of escape velocity and escape energy in the resulting process. They also predicted that 
the motion of charged particle after collision will be chaotic due to the presence of magnetic field 
and strong gravitational field. Although the chances of collision between two particles around 
black hole, in general are feeble. However the process itself is important to describe the ejection of 
particles from the vicinity of black holes. Later on Hussain et al 24] investigated a similar scenario 


for a slowly rotating Kerr black hole and discussed the conditions of escape for the particle. Jamil 
et al 25] investigated a similar scenario for a Schwarzschild black hole surrounded by quintessence. 

Q] obtained both static and non-static black hole solutions in the scalar- 


Recently, Moffat 

tensor-vector modified gravity, the theory which he himself proposed 


scribes several astronomical and cosmological observations such as galaxy rotation curves 
gravitational lensing 28]. The modified gravitational field equations are given by flQ: 

1 


261 ]. The theory fairly de- 
27] and 

r 


R, 


/iv 2 9^R — 8ttG 


(1) 


where 




1 


9fiv 1 


><T0 


4vrV r 

and B^ = (p Ujfl — where is a vector field with the source charge Q = \ZoGnM (see details 
below). The role of this vector field is to produce large scale repulsive gravity which can cause 
accelerated cosmic expansion. Further the vacuum field equations are 


B% = 0 , 


B 




= 0 


( 2 ) 


where ; denotes the covariant derivative operation. 

To solve the above system of equations, an ansatz for the static and spherically symmetric 
solution is assumed of the form 

ds 2 = g(r)dt 2 — g(r) _1 dr 2 — r 2 (d0 2 + sin 2 6d<j) 2 ). (3) 

The calculation yields 


^ , n GM G n GM 2 

g{r) = 1 - 2-b a -5-, 


G 


G — (jr7\r( 1 T o) , Q — kAL, k — i \JoGj\f , oc — —— — 1, 

Gn 


( 4 ) 













and therefore we have 


Q = ± \J olGj^M. 


( 5 ) 


Note that a is a free parameter of the theory, hence it yields a variable gravitational constant. Here 
M and Q is respectively mass and electric charge of the black hole and Gn is the Newton’s gravi¬ 
tational constant. In metric ([3]), positive value of Q is chosen to maintain repulsive gravitational 
force as it is necessary to describe the stable star. (For a stable star, the gravitational attraction 
should be balanced by the repulsive gravity). This is a static, spherically symmetric, point particle 


31 


, 33]. For a = 0, 


solution of a electrically charged body like Reissner-Nordstrom black hole 
metric ( 02 ) reduces to the Schwarzschild metric (which is also the general relativistic limit). Like 
Kerr 133] and Reissner Nordstbrm metrics it has two horizons: 


r± = GjyM ^1 + a ± (1 + a) 2 ^ . 


( 6 ) 


Eq. ([6]) corresponds to g{r) = 0 and will reduce to Schwarzschild event horizon for a = 0. The 
metric ([3]) is invariant under time translation and rotation around symmetry axis. Thus the Killing 


vectors equations are 


34]: 


= d t > CA = A 


(4>r^ 


( 7 ) 


which will give the constants of motion, where = (1, 0, 0, 0) and = (0, 0, 0, 1). The 
conserved quantities corresponding to these Killing vectors are the total energy (per unit mass) £ 
and azimuthal angular momentum L z (per unit mass) of the moving particle at infinity. Motion 
of a neutral particle moving in MOG background is described by the Lagrangian density j^ ]. 


£ = 7>9iwX ii x 1 '. 


( 8 ) 


From (|4]) and (|5J) we can say that t and 4> are the cyclic coordinates. There exists constants of motion 
corresponding to these cyclic coordinates i.e. total energy and azimuthal angular momentum. We 
have calculated these integral of motion by using the Euler-Lagrange equation 

dC 


d / dC \ 
dr V dx^ ) dx^ 


= 0 . 


( 9 ) 


Therefore, using equation ([9|) for t and 0 we have 
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The over dot denotes the differentiation with respect to proper time r throughout the calculations. 
Considering the planar motion of the particle i.e. for 8 = 7r/2 the normalization condition = 1, 
gives 


t 2 = e 2 - u eS , 


t/eff = (l — 2-h 


GM G n GM 2 


L 2 


a- 


)( 1 + dt)- 


( 12 ) 


The extreme values of the effective potential correspond to = 0. It occurs at r = 6 M for a 


Schwarzschild black 
effective potential 


role [23]. The point where the ISCO exists is the convolution point of the 


36]. In the present case, the ISCO occurs at 


r n = 


aGM 2 G N + L 2 
3 GM 


- V 2 


- 3 GM) + aGM 2 G N ) ( L(3GM + L) + aGM 2 G N ) S j 


( 


3 GM I - 2 a 3 G 3 M 6 G% + 3GL 4 M 2 (9 G - 2 aG N ) - 2L 6 - 3 aG 2 L 2 M A G N (2 aG N + 9 G) 


+3a/3 (g 3 L 2 M 4 (^~ 9GL 6 + 108 G 3 L A M 2 + aG N (8L d - 126 G 2 L A M 
+aGM 2 G N (2AL 4 - 9 G 2 L 2 M 2 + 8 aGM 2 G N ( aGM 2 G N + 3 L 2 ))) 


- 2 a 3 G 3 M 6 G 3 N + 3 GL*M 2 (9 G - 2 aG N ) - 2L b - 3aG 2 L 2 M*G n (2 aG N + 9 G) 


' 4 71 ,r2 


1 \ 3 

2 


2 t 2 71 ,r4/ 


1 -1 


3^2 GM 

+3\/3 (g 3 L 2 M 4 ^ - 9 GL 6 + 108G 3 L 4 M 2 + aG N (8L e - 126 G 2 L 4 M 


+aGM 2 G N (2AL A - 9G 2 L 2 M 2 + 8 aGM 2 G N ( aGM 2 G N + 3L 2 ))) 


(13) 


The critical energy and the azimuthal angular momentum of the particle corresponding to ISCO 


are 


£0 = 


(aGM 2 G N + r(r - 2GM)) 2 
r 2 ( 2aGM 2 G N + r(r - 3 GM)) ’ 


r — 

*-' 7,0 - 


y/GMr 3 - aGM 2 r 2 G N 


(14) 


(15) 


s/2aGM 2 G N - 3 GMr + r 2 ' 

We consider the case that an incoming particle collides with the orbiting particle at the ISCO, 
so that after collision it will move within a new plane tilted with respect to previous equatorial 
plane. However, to study the dynamics of a particle it is convenient to use the fact that if the 
initial position and the tangent vector of the trajectory of the particle lies on a plane that contain 
the center of the body, then the entire trajectory lies on this plane. After collision, there are 
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three possible cases depending on the collision mechanism : (i) bound motion (ii) capture by BH 
(iii) escape to infinity. For a little change in energy and angular momentum, orbit of the particle 
alters very slightly. While for large changes it may escape to infinity or capture by BH depending 
upon the nature of change in path. After collision, the particle will no longer remain in the same 
equatorial plane, so further discussion would be dealt with respect to new plane. But note that 
due to spherical symmetry all equatorial planes are equivalent. Due to collision particle should 
have new constants of motion £. L 2 and L z . For simplification of our problem we consider the 
case of collision when (i) the azimuthal angular momentum is invariant during collision ( ii ) initial 
radial velocity does not change. These conditions imply that only energy of the particle will change 
hence its motion would be determined by considering only the change in the energy. After collision 
particle acquires an escape velocity v, in orthogonal direction of the equatorial plane as explained 
in [jjlj]. The angular momentum and energy of the particle after collision becomes (at the equatorial 


plan 0 = 


T2 


L 2 = r 2 v 2 + ~ zo 


j 2 2 2 i t 2 

~ 2 n i ^ 1 cT T L zo 

sm 0 


(16) 


Here v = —r o 0 o and 0 o is the initial polar angular velocity of the particle and r 0 is the radius of 


ISCO. This velocity should be in orthogonal direction to the equatorial plane 


37]. Further 


f 2GM +f OT !) „ 2 + c 


(17) 


where £ 0 is defined in Eq. (|14|> . It is clear that these values of angular momentum (11611 and energy 
(] 171) are larger than their values (|14l) and ()15l) before collision. Also Eq. (|17l) shows that as r -> oo, 
£ —>• £ a —>• 1. So for £ = £ > 1 particle will have unbound motion. In other words the particle can 
not escape to infinity if £ < 1. Hence the necessary condition for a particle to escape to infinity 
after collision is £ > 1 or 


v > 


y/r (2 GM (L 2 + r 2 ) - L 2 r) - aGM 2 G N (L 2 + r 2 ) 


y/r 2 ( aGM 2 G N + r(r - 2 GM)) 

The last expression for velocity v is obtained by solving equation (fT71) at £ new = 1. 


(18) 


III. MOTION OF A CHARGED PARTICLE AROUND A WEAKLY MAGNETIZED 

SMOG BLACK HOLE 

A. Magnetic Field in the Vicinity of BH 


le mag netic coupling (MC) process is responsible for attraction of black hole with its accretion 


T 

disc [171.119 


384401]. According to this process angular momentum and energy are transferred from 
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a black hole to its surrounding disc. The process of MC provides the relation between the strength 
of magnetic field at the black hole horizon and its mass M and the rate of accretion M 
relation is as follows: 


41]. This 


Bh = —\/2 m v Mc. 
fh v 


(19) 


Here black hole horizon r/, is given by ([6]) and m p is the magnetization parameter which indicates 
the relative power of process of MC with respect to disc accretion. If disc accretion dominant over 
MC process then m p < 1 and if MC process is dominant over disc accretion then m p > 1 while 
m p = 1 correspond to equipartition of these two processes. The magnetic field expression is given 


by 


m 


B h = 




x 10 


7 . 35+.450 


( 20 ) 


Magnetic field strength at the horizon of MOG-BH is 4.93 x 10 8 Gauss for m p = 1, a = 0.1, 6 = 0.5 
and Vb = 300. 

Similar effects of particle collision with high center of mass energy in the vicinity of black hole 
can also possible if the black hole is non rotating provided there exists a magnetic field in its 
surrounding. Since, there exists theoretical and experimental evidence that such a magnetic field 
should exist in the surrounding of black hole 2l|]-[29|]. We assume that this magnetic field is weak 


and its energy and angular momentum does not effect the background geometry of the black hole. 
The above mentioned condition satisfy for a black hole of mass M if magnetic field strength holds 
the condition 


B < B — 


1 




G 


M, 


O 


M 


10 


19 M 0 

M 


( 21 ) 


These kind of black holes are called weakly magnetized. 


B. Magnetic Field Calculation 


In this section we explore how does the presence of a magnetic field in the BH exterior stimulate 
the motion of a charged particle. Firstly, we will calculate the magnetic field in the vicinity of black 
hole by following the procedure as given by { 23 , 24] 


The Killing vector equation is 


42] 


= 0 , 


( 22 ) 
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where is a Killing vector. Eq. (1221) corresponds to the Maxwell equation for 4-potential A^ in 
the Lorentz gauge A M ;/i = 0. Defining A M as 43, 44] 

'aGM „ „ 2T 


A» = 


-.0,0, f). 


(23) 


Here magnetic field is defined as 


m 


B>‘ = --e^F Xa 


(24) 


where 


livXcr 


e ° 123 — 9 — det(g^ v ). 


(25) 


Here is the Levi Civita symbol, and F. w is the Maxwell tensor, defined as 


F — A —A —A —A 

-L /_Ll/ - //. Ji-u jy - Jljj-ii. -FL 


t 








For a local observer at rest, the component of four velocity are 

«'* = («*, 0,0,0) = (<?(r)^,0,0,o). 


(26) 


(27) 


Here we assume u l > 0, to signify the ‘forward-in-time’ condition. By equations 
magnetic field 4-vector becomes 


B» 


0 ,Bg(r 


1 

I 2 cos ( 


Bg(r 


■ sin 6,0 


H23H2ZD the 

(28) 


The ma gne tic field is along with the z-axis at spatial infinity, and we assume that it is directed 


upward 


451 ]. At equatorial plane only B e component will survive. In figure [T] we have plotted B 6 


as a function of r. It is decreasing initially and then become almost constant for large r (away 
from the black hole). So, it is homogeneous at r —>• oo. 



FIG. 1: In this figure we have plotted the magnetic field B e vs r for a = 0.2. 
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C. Dynamical Equations 


The Lagrangian of the moving particle of rest mass m and electric charge q in the vicinity of 
MOG-BH is 


£ = + —A^. 

2 m 


Using Euler-Lagrange equation (1271) for t and cj) we get 

. , / ■ eaGM \ 

s(r) (‘ +— r - )“ ’ 


Here, e = is the specific charge of a particle and 


L, 


r 2 sin 2 0 


+ B , 


where 


B = 


eB 


Using the normalization condition = 1, we obtain 

1 = g(r)t 2 -—r 2 — 7,2 ~ 7,2 S7n2 

9{r) 

By using equations (1301) and (13TT) in (l33l) and choosing 0 = f we have 

r 2 + U eS = £ 2 , 

and then effective potential is 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


£± = U eS ± = ± \ g(r) [l + r 2 (^£ - B) 2 


According to equation (fl6l) . after collision L z —>• L. Hence the effective potential reduces to 

eaGM 


(35) 


£± = UeS± = 


± \lg(r) 1 + r 2 (^-B) 2 


(36) 


Putting the value of L from equation (1161) and £± = 1 in equation (1361) and then solving for v we 
get 

1 


v = 


GMr 2 ( Br 2 - L) ( aMG N - 2 r) + Br 6 - Lr 4 


r 4 ( aGM 2 G N + r(r - 2 GM)) 

±sjGMr & (aGM 2 G N + r(r - 2GM)) ( a 2 GMe 2 - aMG N + r(2 - 2ae)) 


(37) 
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For a particle to escape from the black hole’s vicinity, its velocity should be greater or equal to v. 
A charged particle moving in an external electromagnetic field F /IV obeys the equation of motion: 


+ r»x v x a = —F£x° 
m 


(38) 


The dynamical equations for 9 and r respectively are 


;; — 2 • LlcOsO 2 . „ 

8 = — r6 H— f -—o- B z sm 9 cos ( 

r r 4 sin 4 9 


r = 7 — 9 2 (yr 2 — rg(r )) — 


L 2 


r 2 sin 2 9 


(7- 




—B 2 sin 2 9^ — 3rg(r ) + r 2 7 ^ + BL ^27 — —^, 


(39) 


(40) 


where 7 = . For 9 = | the equation (l39l) is satisfied obviously and equation (1401) becomes 

L 2 


r = ry — 


7 - £(7 - - B 2 (—3rg(r) + r 2 7 ) + Bl(2 7 - . 


4g(r) 
r 


(41) 


We have solved equation (1411) numerically by using the built in command of Mathematica NDSolve 
for a = 0.2, B = 0.3 and L z = 2 and plotted the solution in figure [2j In this figure the upper curve 
represents r(r), middle is for r (radial velocity) and the lower one is for r. r is the radial velocity 
which represents the escape behavior as it is increasing with the increase of r, and r is the radial 
acceleration which is analogous to the force exerted on the particle in the radial direction after the 
collision. 


r(T)r'(r)r"(T) 



FIG. 2: Interpolating function r(r) as the solution of Eq. (l4lT) for L z = 2, B = 0.3, and a = 0.2. Here the 
large bold dashed curve represents rr, solid curve is for r{r) and Short dashed curve correspond to r(r). 
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IV. FORCE ON A CHARGED PARTICLE IN THE VICINITY OF MOG-BH 


As we have already calculated the effective potential for MOG-BH, we can also compute the 


effective force on a neutral and a charged particle by 

p _ 1 dU e s 

2 dr 


(42) 


p _ -G n M (3 L 2 + r 2 ) + L 2 r + aG N M ((a + 1 )G N M (2 L 2 + r 2 ) - r (3 L 2 + r 2 )) 

It can be seen from equation (1431) that the force due to scalar tensor-vector gravity is repulsive if 
(a + 1 )G n M (2 L 2 + r 2 ) > -r (3L 2 + r 2 ). 


r _ L 2 r - GM (3 L 2 + r 2 ) | Q 2 (2L 2 + r 2 ) 


(44) 


Equation (1441) represents the effective force for RN-BH and the force due to charge of BH is repulsive 
with out any condition. 

We are studying the dynamics of a neutral and a charged particle in the surrounding of MOG- 
BH where the scalar tensor-vector field produces a repulsive gravitational force which prevents a 
particle to fall into singularity [26j]. In figure [3] we are comparing the effective force on a particle 


in the vicinity of MOG-BH with Schwarzschild black hole. We deduce from this figure that the 
repulsion to reach the singularity is more for a = 0.2 as compared to a = 0. 



r 


FIG. 3: Here we have plotted effective force as function of r for different values of a. 
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To study the behavior of force against magnetic field B we have plotted the force F as a function 
of r for different value of B in figure [4l It can be seen from this figure, due to large magnetic field 
strength B, force on the particle is more attractive as compared to small B. 



FIG. 4: Here we have plotted effective force against r for different values of magnetic field B. 


V. COMPARISON OF GEODESICS IN THE VICINITY OF MOG-BH VS SBH 

Geodesics of a Neutral Particle Moving Around a Schwarzschild BH: 

Geodesics of a particle approaching toward or away from BH could be obtained by using Eqs. cnD 
and (fl2l) together. We have: 

dt £ 1 , . 

~r~ — -f- 7 r 7 , (45) 

dr g{r) ^£2 _ UeS 

Here U e g correspond to a neutral particle given by (fT2l) . where positive root gives the path of the 
particle out-going from the BH, and negative root gives the path of an ingoing particle. Let us 
consider the particle which is coming from infinity, initially at rest, approaches the BH, so setting 
£ = 1, B = 0, a = 0 and L z = 2 in Eq. (PT5l) we plot the geodesics in (r, t) coordinates, see figure 

0 

Geodesics of a Charged Particle Moving Around a MOG-BH: 

Geodesics of the particles approaching the MOG-BH could be obtained by using Eqs. (1301) and 
(1341) together we obtain: 

j £g(r )-1 - 

dr V £ 2 “ U eS 


(46) 
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0 5 10 15 20 

distanc - r 


FIG. 5: Here we have plotted the geodesic equations (l45l) against r for £ = 1, L z = 2, M = 1. 

Here J7 e ff corresponds to a charged particle given by equation (l36l) where positive and negative signs 
give the path of the outgoing and ingoing particle respectively. Setting 5 = 1, L z = 2, M = 1, in 
Eq. (1461) we get the geodesics which are bounded by the boundaries r = r c and the outer horizon of 
the BH, plotted in hgure[6l In this figure, for thick curve we consider parameter a = 0.1, magnetic 
field strength B = 0.1, and charge of a particle e = 0.5 and thin curve correspond to e = 1, B = 0.5 
and a = 0.5. Figure [6] shows that if the strength of magnetic held B is higher, then charged particle 
can reach arbitrary close to black hole as compared to smaller value of B. 



0 2 4 6 8 

distanc - r 


FIG. 6: Here we have plotted the geodesic equations (l46ll against r for £ = 1, L z = 2, M = 1. 
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VI. STABILITY OF ORBITS 


Lyapunov exponents are the measurements of the rate of convergence or divergence of the nearby 
trajectories in the phase space. It is highly sensitive to the initial conditions. Its positive value is 
the indication of divergence among the nearby trajectories. Therefore we can check the stability 
of orbits by the Lyapunov exponent A 47[. It is given by 


A = 



(47) 


where r a is the ISCO of the particle moving around BH. 
instability is defined as 


In 4 


a critical component for the 




(48) 


Here T\ is the instability time scale. Time period for a circular orbits can be obtained from equation 
m in the absence of magnetic held as 


T = 


^r 2 0 

L z 


(49) 


and 


in the presence of magnetic held, form equation (1311) we have 


T =1 


2tt ri 


-T 4- Rr2 


(50) 


Here r 0 is the radius of circular orbit. We have plotted Lyapunov exponent as a function of 
magnetic held B in figure [71 It shows that the larger the magnetic held strength the lesser will be 
A. From figure [7] and equation (1501) one can say that the instability of circular orbits is more for a 
Schwarzschild black hole in comparison with the black hole immersed in a magnetic held. In hgure 
[8] we have also plotted A against a it decreases by the increase of a. Hence, the stability of circular 
orbits is less for the Schwarzschild black hole as compared to the black hole with non zero a. 
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FIG. 7: In this figure we have plotted the Lyapunov exponent as a function of magnetic field B for L = 6, 
a = 0.5, £ = 1, and e = 1. 



FIG. 8: Plotting of the Lyapunov exponent as a function of parameter a. 


VII. BEHAVIOR OF EFFECTIVE POTENTIAL 


Effective potential largely depends on goo component of metric. Therefore, before discussing 
the behaviour of effective potential we will compare the goo component of MOG-BH metric with 
RN-BH metric as it look similar. For RN-BH metric, 


1 _ 2 M + Q 


9oo = 


(51) 
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for MOG-BH from equation Q we have 

( 2 M , Q 2 \ , s 

9oo = ( 1--- 2y/a-^ + y/a{l + a)-^-J. (52) 

One can see that equation (1511) contain only Q 2 term while equation (1521) contain Q and Q 2 . This 
difference leads to large change between the behaviour of effective potential of these two black 
holes. The goo component of MOG-BH metric also contain a parameter a which will create a main 
difference between the behaviour of potential. 

In this section we plot the effective potential and graphically explain the conditions on the 
energy of the particle required for escape to infinity or for bound motion around MOG-BH. In 
figure [9] and figure [12] we have plotted t/_eff and U .|_eff respectively correspond to equation (l36|) . 
Here we will discuss the long distance and short distance behaviour of equation (1361) . One can 
see that for small value of r the square root term will dominate but for large r the term which is 
proportional to ^ will dominate. For U+ we get some maximum value which we have mentioned in 
figure M as U max . For a particle to fall into the black hole its energy should be greater then U max 
otherwise it will bounce back to infinity or to some stable orbit. In case of U-, particle can never 
cross the barrier as shown in the figure [9j Hence it cannot fall into the singularity. 



r 

FIG. 9: In this figure we have plotted the effective potential U+ against radial coordinate r. 

In figure fill different regions of effective potential which correspond to escape and bound motion 
of the particle are shown. In figure [TTl /3 corresponds to region of stable orbits. If the particle has 
energy equal or greater then 6, then it will fall into the black hole. U max and U m i n correspond to 
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FIG. 10: In this figure we have plotted the effective potential U+ as a function of r. 


unstable and ISCO orbits respectively. If the energy of the the particle is equal or less than (3 then 
it will reside in one of the stable orbits. Particle having energy greater than /3 but less or equal 
to 5 can have two possibilities: either it will escape to infinity or captured by the black hole. The 
local minima might be related to the change in behavior of the effective potential corresponds to 


BH (black brane transition) as observed in 48]]. 


In figure IT?1 we are comparing the effective potential of Schwarzschild-BH, RN-BH and MOG- 
BH. It can be seen from this figure that for large r all the potentials approach to 1. Hence it can 
be concluded that the minimum energy required for a particle to escape from black hole vicinity 
is £ = 1 as we have mentioned before. Further the maxima for the effective potential of RN-BH is 
greater in comparison with the maxima of effective potential of Schwarzschild-BH and MOG-BH. 
Particle will be captured if it has energy greater than this maxima otherwise it will move back 
to infinity or may be reside in some stable orbit. Therefore, we can say that the possibility for 
a particle to escape from the vicinity of a black hole or to reside in some stable orbit is high in 
case of MOG-BH with a = 1.8 in comparison to RN-BH and MOG-BH with a = 0.8. Particle can 
reach to the singularity depending upon its energy, but it cannot reach the singularity for a = 1.8 
as shown by the curve, in figure fl2l Note that a > 1, corresponds to a naked singularity. 

In figure [13] we are comparing the effective potentials in the presence with the absence of 
magnetic field. In this figure u ma xi and u ma x 2 correspond to unstable orbits while u m i n \ and u m i n 2 
is refers to ISCOs. One can notice that in the presence of magnetic field, minima of the effective 
potential is shifted towards the horizon and width of stable region is also increased in comparison 
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r 


FIG. 11: In this figure different regions of effective potential which correspond to escape and bound motion 
of the particle are shown. Here /3 correspond to stable orbits for b = 0.5. 



r 


FIG. 12: Comparison of effective potential U e s as a function of r. In this figure we are comparing the 
effective potential of RN-BH and MOG-BH. 


with the case when the magnetic is absent. This is in agreement with [24, |4J]. Therefore we can 




























say that magnetic field act as to increase the stability of the orbits. 



r 


FIG. 13: Behavior of effective potentials with and without magnetic field versus r (a comparison). 

In figure [14] we have plotted the effective potential as a function of r for different value of angular 
momentum L z . One can see that for large value of L z . the effective potential has local minima 
and maxima which correspond to stable and unstable circular orbits respectively. Hence we can 
say that the particle with greater value of L z would have more possibility to reside in the stable 
orbits in comparison with lesser value of L z . 

To study the behaviour of effective potential against e, we have plotted effective potential as a 
function of radial coordinate r for different values of e in figure [15j It can be seen that larger the 
value of e greater will be the maxima of potential. Therefore, large value of e would correspond 
to easy escape and vice versa. It can be seen from the detailed analysis of the effective potential 
and the black hole physics that the scalar tensor-vector modified gravity differs from the Einstein’s 
theory of gravity a lot at shorter distance and it becomes similar at long distance. 

VIII. TRAJECTORIES OF ESCAPE VELOCITY 


For the angular variable we have, 
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Effective Potential For Different value of L 



r 


FIG. 14: Comparison of effective potentials with respect to L z as a function of r. 



r 


FIG. 15: Here we have plotted effective potential as function of r for different values of e. 


f the left hand side of equation (1531) is negative then the Lorentz force on the particle is attractive 
4s| . The motion of the charged particle is in clockwise direction. Lorentz force is repulsive if left 
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land side of (l53l) is positive. We are not going to the detail here because it is already discussed in 


36, 


49]. Our concern is only about the action of magnetic field on the charged particle. Therefore, 


magnetic field may deform the oscillatory motion, so, the greater the strength of magnetic field, 
the larger will be the deformation of the orbit. Hence we can conclude that larger the strength of 
magnetic field, easy for a particle to escape from the ISCO. 

We have also plotted in figure fTHl escape velocity as a function of radial coordinate for different 
values of magnetic field parameter b. It can be seen that the escape velocity of the particle increases 
as the magnetic field strength increases but it becomes almost constant just like magnetic field, 
away from the BH . As the magnetic field is strong near the BH, therefore, we can conclude that 
presence of magnetic field will provide more energy to partic e, so that it might easily escape from 


the vicinity of BH. These conclusions are consistent with 


HQ. 



FIG. 16: Escape velocity ( v) as a function of r for different values of magnetic field B. 

Figure [TT] represents the escape velocity against radius r for different values of angular momen¬ 
tum L z . From figure [171 we can say that the possibility for a particle to escape having large angular 
momentum is small. We have plotted escape velocity for different values of a in figure [18] and [l8l 
a = 0, correspond to S-BH (Schwarzschild Black Hole) and a = 1 correspond to RN-BH. 

In figure flHI we are comparing the escape velocity of a particle moving around the S-BH, RN- 
BH and MOG-BH. Note that the difference between the velocities is larger near the black hole 
(initially) and it becomes almost same away from the black hole. Therefore, we can conclude that 
the effect of the charge of black hole on the motion of the particle is large while it reduces as 
particle moves away from it. One can see that for large r escape velocity is same for all values of 
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FIG. 17: Escape velocity (v) against r for different values of angular momentum L 
a but for small r lesser value of a correspond to greater value of escape velocity and vice versa. 


Escape Velocity for Different Values of a 



FIG. 18: Escape velocity ( v ) against r for different values of parameter a. 
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IX. SUMMARY AND CONCLUSION 

We have investigated and compare the dynamics of a charged and a neutral particle in the 
vicinity of S-BH, RN-BH and MOG-BH. Geodesics of a neutral and a charged particle in the 
vicinity of a MOG-BH are shown in figures [5] and [6l We see that for a charged particle there are 
two boundaries on the geodesics, r = and r = r c unlike the formal case in which a neutral 
particle comes from infinity and goes back to infinity before reaching the horizon, r = of the 
BH. We discussed the effective potential behavior in details regarding the stability of the orbits 
of the particle. We further discussed the energy condition for the particle, when it will escape 
or its motion remains bound. Expressions for the escape velocity of the particle moving around 
MOG-BH and for magnetic field, present in the vicinity of BH due to plasma, is derived in this 
work. More importantly a comparison is done for effective potentials, obtained in the presence and 
absence of magnetic field. It is found that presence of magnetic field enhance the stability of the 
orbits of the moving particles, due to the presence of it width of stable region in contrast to which 
we obtained in the case when magnetic field is absent. We have also done the comparison of the 
effective potentials among the RN-BH, S-BH and MOG-BH. Further we studied the stability by 
Lyapunov exponent against magnetic field and a parameter a. We conclude that stability of orbits 
would increase due to the presence of vector field considered in MOG. We deduce that the particle 
have to face more repulsion to reach the singularity due to the presence of vector field as considered 
in MOG. But the presence of magnetic field might increase the attractive force. It is found that a 
particle with large value of angular momentum L z would have greater possibility to reside in the 
stable orbits in comparison with a particle with lesser value of it. Therefore, escape velocity is less 
correspond to particle with large value of L z . Effects of magnetic field and parameter a on escape 
velocity are also investigated graphically. It is concluded that presence of magnetic field might 
provide sufficient energy to particle to escape easily from the surrounding of BH. 
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